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Collective modes of bilayered superconducting superlattices (e.g., YBCO) are investigated within 
the conserving gauge-invariant ladder diagram approximation including both the nearest interlayer 
single electron tunneling and the Josephson-type Cooper pair tunneling. By calculating the density- 
density response function including Coulomb and pairing interactions, we examine the two collective 
mode branches corresponding to the in-phase and out-of-phase charge fluctuations between the two 
layers in the unit cell. The out-of-phase collective mode develops a long wavelength plasmon gap 
whose magnitude depends on the tunneling strength with the mode dispersions being insensitive 
to the specific tunneling mechanism (i.e., single electron or Josephson). We also show that in 
the presence of tunneling the oscillator strength of the out-of-phase mode overwhelms that of the 
in-phase-mode at fcy = and finite k^, where k^ and fcy are respectively the mode wave vectors 
perpendicular and along the layer. We discuss the possible experimental observability of the phase 
fluctuation modes in the context of our theoretical results for the mode dispersion and spectral 
weight. 
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I. INTRODUCTION 

In contrast to bulk isotropic superconductors, where 
longitudinal collective modes (i.e., plasmons) associated 
with the density fluctuation response is virtually of no 
particular interest or significance in the context of super- 
conducting properties, there has been substantial recent 
theoretical interest in the longitudinal collec.tive mode 
spectra of layered high- Tc superconductoraJlj. This in- 
terest arises primarily from the highly anisotropic two 
dimensional layered structure of these materials which, 
in principle, allow for sub-gap plasmon modes resid- 
ing inside the superconducting gap in the low wave 
vector regime. ^ This gives rise to interesting collective 
mode behaviorEiLj in layered anisotropic superconduc- 
tors which have no analogs in bulk isotropic supercon- 
ductors. In this paper we consider the effect of having 
multilayer complex unit cells, as existing in YBCO and 
BISCO high-Tc superconductor materials, on the lon- 
gitudinal electronic collective mode spectrum. We find 
a number of collective modes arising from the complex 
unit cell structure, and comment on their possible ex- 
perimental relevance. One of our goals is to critically 
assess whether observable electronic collective mode be- 
havior could shed some light on the interesting and un- 
usual mechanism producing high Tc superconductivity in 
these materials. The other goal is to predict novel collec- 
tive mode behavior peculiar to layered superconductors 
with no analogs in bulk systems. 

The collective mode spectrum is characterized by the 
energy dispersion {h = 1 throughout this paper) rela- 



tion Lv = w(fc||,/cz), which we calculate in this paper, 
where fcy = |k||| is the two dimensional wave vector 
in the so-called a-b plane (along the layer) and is 
the wave vector along the c-axis, then kz — |k|cos^, 
fc|| = |k|sin0. Because of the strong a-b plane versus 
c-axis anisotropy in these materials, the dependence of 
the collective mode frequency on fcy and k^ is very dif- 
ferent. [We ignore any anisotropy, which is invariably 
rather small, in the a-b plane and assume intralayer pla- 
nar isotropy, i.e., LoCk^^jk^) = uj(k\\, kz).] The structural 
model we employ considers the layered superconductor 
to be a one dimensional superlattice along the z direc- 
tion (c-axis) composed of a periodic system of bilayer 
unit cells with an intracell layer separation of c and a 
superlattice period of d (> c). The two active layers sep- 
arated by a distance c within each unit cell are taken to 
be identical and are assumed to be planar two dimen- 
sional electron gas (2D EG) systems of charge density 
Us per unit area and zero layer thickness each. In most 
of our calculations presented in this paper the intercell 
electron hopping (or tunneling) between neighboring unit 
cells (separated by a distance d) is neglected (i.e., we ne- 
glect any superlattice band width along the z direction), 
but we critically examine the effect of intracell electron 
hopping between the two layers within each unit cell on 
the collective mode dispersion. We comment upon the 
effect of a finite mtercell hopping in the conclusion of 
this article. We include in our theory the long range (in- 
tracell and intercell) Coulomb interaction among all the 
layers. This long range Coulomb interaction, which cou- 
ples all the layers, is of great importance in determining 
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the collective mode spectrum. We also include in our 
theory of collective mode dispersion the effect of the su- 
perconducting pairing interaction, assumed in our model 
to be a short-range in-plane attractive interaction of the 
BCS-Fermi liquid type, which is treated in a fully gauge 
invariant Nambu-Gorkov formalisKU, Our work is thus a 
generalization of the earlier workldia by Fertig and Das 
Sarma, and by Hwang and Das Sarma (who considered 
only the monolayer superconducting superlattice situa- 
tion with only a single layer per unit cell) to a complex 
unit cell situation with two layers per unit cell. To keep 
the situation simple we will consider only the s-wave gap 
symmetrytl, which, according to ref. |^ gives a very good 
account of the collective mode dispersion even for the 
d-wave case except at very large wavCr-vectors. Follow- 
ing the work of Fertig and Das Sarmald there has bf fQi 
a great deal of theoretical and experimental workcTtHl 
on the electronic collective mode properties in layered 
superconducting materials, but the specific issue consid- 
ered in this paper has not earlier been discussed in the 
literature for a multilayer superconducting system. It 
should also be pointed out that, while the focus of our 
work is the collective mode behavior in layered high-Tc 
cuprate superconductors (which are intrinsic superlattice 
systems due to their highly anisotropic crystal structure 
with CuO layers), our results equally well describe ar- 
tificial superconducting superlattices made of multilayer 
metallic structures provided fc|| and kz are good wave 
vectors in the system. 

The collective mode dispersion in bilayered supercon- 
ducting superlattices is quite complicated. There are es- 
sentially two different branches of long wavelength collec- 
tive modes: in-phase (w+) modes and out-of-phase (w-) 
modes, depending on whether the electron density fluc- 
tuations in the two layers are in-phase or out-of-phase. 
Each of these collective modes disperses as a function of 
wave vector, showing strong anisotropy in fcy and kz dis- 
persion. In particular, the limits {kz — 0, fcy — > 0) and 
(kz — !■ 0, fc|| — 0) are not equivalent because the kz = 
three dimensional limit is singular. For kz = the in- 
phase collective mode is a gapped three dimensional 
plasma mode at long wavelengths {kz — 0, fcy — > 0) by 
virtue of the Higgs mechanism arising from the long range 
Coulomb interaction coupling all the layers. This mode 
characterizes the long wavelength in-phase charge fluc- 
tuations of all the layers. For non-zero kz the w+ mode 
vanishes at long wavelengths (fc|| 0) because at finite 
kz the system is essentially two dimensional. The out-of- 
phase LU- collective mode branch arises purely from the 
bilayer character of the system and indicates the out-of- 
phase density fluctuations in the two layers. In the ab- 
sence of any interlayer hopping (either intracell and inter- 
cell) the uj- mode is purely acoustic in nature vanishing 
at long wavelengths (/cy 0) as uj-{kz,ki\ ^ 0) ~ 0(A:j| ) 
independent of the value of kz ■ For finite interlayer tun- 
neling oj- exhibits a tunneling gap at fcy = 0. The Higgs 
gap for uj+{kz = 0,k^\ —f 0) is not qualitatively affected 
by intracell interlayer tunneling because the three dimen- 



sional plasma energy is usually substantially larger then 
the tunneling energy. 

Note that, in the absence of any intracell and inter- 
cell tunneling, both in-phase and out-of-phase collective 
mode branches lie below the superconducting energy gap 
for small [except for the LU-i-{kz — 0) mode which is 
pushed up to the three dimensional plasma frequency]. 
This remains true even for weak intracell and intercell 
tunnelings, and in this paper we concentrate mainly on 
this long wavelength "below gap" regime where the phase 
fluctuation modes could possibly lie in the superconduct- 
ing gap. For simplicity we also restrict ourselves to s- 
wave gap symmetry of the superconducting order param- 
eter. This approximation may at first sight appear to be 
unusually restrictive as it seems to rule out the applica- 
bility of our theory to bilayer high- Tc materials (such as 
YBCO, BISCO) which are now widely accepted to have 
rf-wave ground state symmetry. This, however, is not 
the case because at long wavelengths (small fcy), which is 
what we mostly concentrate on, the collective mode spec.- 
trum is insensitive to the order parameter symmetry^, 
and therefore our results apply equally well to high- Tc 
bilayer materials. The modes we predict and their dis- 
persion should most easily be observable via the resonant 
inelastic light scattering spectroscopy, but may also be 
studied via frequency domain far infrared spectroscopy 
using a grating coupler. 

II. THEORY, APPROXIMATIONS, AND 
RESULTS 

In our calculation we assume that the two layers in each 
unit cell can be considered to be 2D EG, and all layers are 
identical, having the same 2D charge density per unit 
area. Two identical layers separated by a distance c in 
each unit cell are strongly coupled through the interlayer 
intracell electron tunneling. The interlayer tunneling is 
between the well-defined CuO layers in high Tc materials. 
The intercell tunneling between different unit cells sepa- 
rated by a distance d (in our model d > c) is neglected at 
first (see section HI for the effect of intercell tunneling) . 
Although we neglect the electron tunneling between dif- 
ferent unit cells, the electrons in all layers are coupled 
via the intercell long range Coulomb potential which we 
keep in our theory. Since the long wavelength plasnia, 
modes are independent of the gap function symmetrya, 
we work in the BCS approximation with s-wave pait 
ing for simplicity. Then, in the Nambu representation^ 
the effective Hamiltonian of a bilayered superconductor 
with 2D quasiparticle energy e{k), a tight-binding coher- 
ent single-electron intracell hopping t{k), and an addi- 
tional coherent intracell Josephson coupling Tj between 
two nearest layers is given by 

H = Ho-iiN + Hi^t + HT,, (1) 

with 
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n,i k 

+EE^w*U^3*k,™, (2) 

n,i k 



-^int — 2 / ^ / ^ Pq,mz^^z,nj(q)P-q,nj; (3) 
ni,mj q 

= E E ^'^ (*k+q,„.T3*k.™) (*k'-q.,»:T3*k',™) , 
n,'i k,k',q 

(4) 

where n, m are the unit cell indices and i, j — 1,2 label 
the two layers within a given unit cell (i = 3 — i). Here, 
^k,Tii and ^'1^ are the column and row vectors 

*k,™ = J^'J";^^ ^ : *k,m = (4,m,T' C-k,ma) > 

(5) 

where cj^. „j ^ (ck,ni,o-) creates (annihilates) an electron 
with wave vector k and spin a in the i-th layer within 
the nth unit cell, and Pq^ni denotes the density operator 
defined by 

Pq,m = E *I+q,m^3^'k,m, (6) 
k 

where £k = /2m — /i (/i being the chemical potential), 
and Ti (j=l,2,3) are the Pauli matrices. In Eq. (||), 
the effective interaction Vni^mjil) contains the long range 
Coulomb interaction coupling all the layers and the short 
range attractive intralayer pairing interaction (giving rise 
to superconductivity in the problem) 

Vni.mjiq) = Vc((7||) exp[-(7|||z„i - Zmj\] + VoS„^rnSt,j (7) 

where Vc((?||) — 27re^/(Kg||) is the two dimensional 
Coulomb interaction and k is the background dielectric 
constant of the system. Vq represents a weak, short- 
ranged attractive intra-layer pairing interaction which 
produces superconductivity, and is a model parameter 
in our theory. 

We should comment on one unusual feature of our 
Hamiltonian defined in Eqs. (|^)-(^). This is the ex- 
istence of both a coherent single-particle hopping term, 
defined by the single-particle hopping amplitude t(k) in 
Eq. (H), and a coherent Cooper pair Josephson tunnel- 
ing term, defined by Tj in Eq. (^. Usually the existence 
of a single-particle hopping t automatically generates an 
equivalent Josephson coupling Tj in the superconduct- 
ing system, and keeping both of them as we do, namely, 
t in the single particle Hamiltonian Hq [Eq. ( p[)] and 
Tj in the two-particle Josephson coupling [Eq. T4|)], is 



redundant. We do, however, wish to investigate sep- 
arately effects of both coherent single particle hopping 
and Josephson coupling along the c-axis on Lhje collective 
mode spectra because of recent suggestionalj of a novel 
interlayer tunneling mechanism for superconductivity in 
cuprates which explicitly postulates i = (in the normal 
state) and Tj ^ (in the superconducting state). Our 
model therefore uncritically includes both t and Tj as dis- 
tinct contributions, and one could think of the interlayer 
Josephson coupling Tj in our model Hamiltonian arising 
from some interlayer pairing interaction not included in 
our model pairing interaction Vq which is exclusively in- 
tralayer in nature. In the following we take t and Tj to be 
independent parameters of our model without worrying 
about their microscopic origins. 

The collective modes of the system are given by the 
poles of the reducible density response function x(k, w). 
We apply the conserving gauge invariant ladder diagram 
approximatioiu'Lj in calculating the density response of 
the system including the effect of the pairing interaction 
induced vertex and self-energy corrections. The density 
response function is defined as 

/>oo 

X(k, c.) - -z / dte^"* < [p(k, i), p(~k, 0)] >, (8) 
Jo 

where p(k, t) is the three dimensional Fourier transform 
of the density operator in the Heisenberg representation. 
Here, k = (A:||, fcz) is the 3D wave vector, where fcz mea- 
sures the wave vector along the z-axis (i.e., the c-direction 
) and is the 2D x-y plane (i.e., a-h plane) wave vec- 
tor. The density response may be written in terms of an 
irreducible polarizability H(k, w) as shown in Fig. 1(a). 
Following Anderson's-arguments for the absence of sin- 
gle particle tunneling)13 we first neglect inter-layer single 
electron tunneling effects (i = 0) and only consider the 

FIG. 1. (a) Diagrammatic representation of the dielectric 
response in terms of the irreducible polarizability 11. Here, Vi 
and V2 axe given in Eq. (|l|), and j = 3 - j. (b) Irreducible 
polarizability used in this calculation. 
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Josephson coupling effect. The polarizability 11 is then 
diagonal in the unit cell index and becomes the corre- 
sponding 2D polarizability matrix, n(k, u) = n(/cj| , u) 



n(fc||,^) 

e(k,t^) ' 



(9) 



where n(fc|| , lo) is the irreducible polarizability matrix for 
a single isolated unit cell, 



^ ' V n2l(fc||,U;) 1122 



(10) 



where Hn, 1122 and ni2, 1121 indicate the intra-layer and 
inter-layer irreducible polarizability, respectively. Within 
our approximation, the inter-layer polarizabilities vanish 
when the single-particle tunneling is neglected. We will 
see that the plasma gap of the out-of-phase mode arises 
entirely from the non- vanishing inter-layer irreducible po- 
larizability. In Eq. (|^) the effective dynamical dielectric 
fimction e(k, lo) is given by 



e(k,w) 



i/(k)n(fc||,c^). 



(11) 



where 1 is a 2 x 2 unit matrix and V^(k) is the effective 
interaction which in our simple model is given by 



T/(k) 



Fi(k) T/2(k)\ 
y2(k) Fi(k) ] 



(12) 



where Vi(k) corresponds to the intralayer interaction 
(« = j) and V2(k) the interlayer interaction (i ^ j), which 
arises entirely from the long-range Coulomb coupling in 
our model, and they are given by 



^i(k) 
^2(k) 



^c(fc||)/(k)- 

^c(fc||)g(k). 



(13) 



where /(k) and .g(k), the superlattice form factors which 
modify the 2D Coulomb interaction due to Coulomb cou- 
pling between all the layers in our multilayer superlattice 
system, are given by 



/(k) 



sinh(fc||(i) 



cosh(/i;j|d) — cos{kzd) ' 



5(k) 



sinh[A;|| {d — c)] 



-ikzd 



sinh(fc||c) ^, 



cosh(fc||(i) — cos{kzd) 



(14) 



(15) 



In order to obtain the collective mode spectrum, it is 
necessary to construct a gauge invariant and number- 
conserving approximation for n(k, w). In the conserving 
gauge invariant ladder diagram approximationEiO the ir- 
reducible polarizability obeys the ladder integral equa- 
tion which is shown diagrammatically in Fig. 1(b). It 
may be written in the form 

Ilij{k,uj) ^ -iTrJ ^^^i^r3Gi(pi, wi) 



where Gi{k,uj) is the i-th layer Green's function with 
self-energy corrections (self-consistent Hartree approxi- 
mation in the Coulomb interaction and self-consistent 
Hatree-Fock approximation in the short-range pairing in- 
teraction) and Fij is a vertex function. The vertex part 
satisfies the linear integral equation 

Fy {pi,k,uj) = T:i5,j+i'Y^ J -j^-^Tsdiq^uJi) 



1=1 ' 



X Fy((7, k,uj)Gi{q -k,,u;- u;,)t3 [Vo5i, + TjS^] , (17) 

where I — 3 — I. In order to solve this vertex function, 
we expand Tij in Pauli matrices 



1=0 



iTl- 



(18) 



Since our model assumes two identical 2D layers in the 
unit cell, we have Fn = F22 = Fq and F12 — F21 — Fb. 
By introducing the polarization function 



(Pqduji 
(27r)3 



T3G{q,uj)TiGiq -k,uji- uj)t3 



the vertex function, Eq. (p^), becomes 



la 
lb 



Pla 

Eib 



( Pjb 

[Pla 



(19) 



(20) 



where 7's are column vectors, = (0, 0, 0, 1), and P is 
a 4 X 4 matrix whose components are given by in Eq. 
(|l9|). Then, the polarizability function 11^ becomes 



i,lT3lij,l 



1=0 



(21) 



(=0 



The poles of the vertex function or polarizability 11 
give the collective mode spectra for the neutral supercon- 
ductor (i.e., neglecting the long range Coulomb coupling 
which couples all the layers) . In the long wavelength limit 
we have two collective modes ( "phasons" ) for the neutral 
system 

t^+(fc) = {vokf [1 + {Vo + Tj)iVo/2] , (22) 



u?{k) = uol + vlk'' [1 + iVo(K) - P/)/2] , (23) 

where = vp/^/^ with Ui? as the Fermi velocity, = 
m/iT is the 2D density of states at the Fermi surface, and 
ol = 16TjAV[^o(Vo^-Pj)] is the tunneling gap induced 



xFij(pi, fc,cj)G'i(fc — — o^i), (16) by the finite Josephson coupling 
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FIG. 2. (a) The plasmon mode (a;±) dispersions in the 
presence of Josephson tunneling for the neutral bilayered 
superconducting superlattice as a function of k\\ for fixed 
fczd = TT. Here, x = Tj/Vo indicates the Josephson tunneling 
strength with respect to the intra-layer pairing interaction. 
Inset shows the ratio of the oscillator strength of a)+ to that of 
cij_. (b) The plasmon mode dispersions for the charged 

system. Inset shows the ratio of the oscillator strength of lj_ 
to that of (c) The a'-(k) band in the superlattice for the 
charged system as a function of in-plane wave vector {k\\(S) in 
the presence of the tunneling. Inset shows the lj± band of 
the bilayer superconducting superlattice. We use parameters 
roughly corresponding to YBCO in these figures: the sheet 
density tIs = 10^*cm~^, effective in-plane mass m — 5mo, 
lattice dielectric constant «: = 4, d = 12^4, and c = SA. 

(Tj 7^ 0). The ti;+ mode corresponds to the in-phase 
motion of the order parameter, or, equivalently the 2-D 
Goldstone-Anderson-Bogoliubov phase fluctuation mode 
due to the spontaneously broken continuous gauge sym- 
metry of the superconducting state. The a;_ mode cor- 
responds to the out-of-ptia£e mode first predicted for a 
two-band superconductorllj, which has recently been cal- 
culated within the time-dependent Hartree-Fork-Gor'kov 
(mean-field) approximatiorjj for a two-layer supercon- 
ductor system. In Fig. ^(a) we show the calculated col- 
lective mode dispersion for different Josephson tunneling 
strengths with respect to the intra-layer pairing interac- 
tion, X — Tj/Vq. When the Josephson tunneling is ab- 
sent, a; = 0, the two phason modes lu± are degenerate and 
have identical dispersion (solid line) . But in the presence 
of finite Josephson tunneling between the nearest layers, 
X ^ 0, the out-of phase mode (uj-) develops a plasma gap 
(luq) depending on the tunneling strength. The in-phase 
mode LUjf- is not affected qualitatively by finite Josephson 
tunneling and remains an acoustic Goldstone mode (i.e., 
cl'+ ^ 0{k) for k ^ Q) although the velocity of the acous- 
tic plasmon does depend on Tj (cf. Eq. (|2^)). In Fig. 
^(a) the inset shows the relative oscillator strength of the 
two phason modes, the ratio of the spectral weight of lj_ 
to that of ix'+ . The ratio decreases as tunneling amplitude 
increases. This is due to the approach of the lj_ mode 
to the pair-breaking excitation region (a; w 2A) at large 
tunneling, which causes decay of the cj- mode to single 
particle excitations, and the strength of the mode trans- 
fers to pair-breaking excitations. These results apply to 
any bilayered neutral superconductors (which, of course, 
do not exist in nature because Coulomb interaction is 
always present in real systems). 

By looking for zeros of the dynamical dielectric func- 
tion defined by Eq. (|l^) we find the collective modes 
of the charged superconducting superlattices. Since the 
two layers within the cell are identical we have IIn = 1122 
and ni2 = 1121, which gives rise to distinct in-phase and 
out-of-phase collective charge density fluctuations of the 
charged superconductor. Coupling of the in-phase (out- 
of-phase) mode of the neutral system via the long range 
Coulomb interaction to the charge density fluctuation of 
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the layers gives rise to the in-phase (out-of-phase) collec- 
tive mode of the charged bilayer system. The dielectric 
function is a matrix, and the zeros of the det[e], which 
define the collective mode spectra, are given by 



det[e] = [l-(nn+ni2)(yi + 1^2)] 

x[i-(nn-ni2)(Vi-y2)] 



0. 



(24) 



In the long wavelength limit Eq. (g^) can be analytically 
solved using Eqs. ( p^ ) - (^ij), and we find two distinct 
collective modes corresponding to the relative phase of 
the charge density fluctuations in the two layers within 
each unit cell: 



a;2(k)=c.^^[/(k) + |5(k)|]. 



+0 ' 



(1 + AT/- Ayo)(cjg + i.gfcV2) 
1 - w2(Al/- Ayo)/6 



(25) 



(26) 



where ojp = (innBe^ / KmY^^ is a three dimensional 
plasma frequency with the effective three-dimensional 
electron density of the double-layered supperlattice ns = 
2ns/ d, and fc^ = fc^ 4- fc| with k = {k\\,kz); = 

Na{Vi - V2) and AFo = N^iVo - Tj)/2. In Fig. |(b) 
we show the calculated charge density mode dispersion 
for fixed k^d = tt as a function of k\\d. Tunneling has 
little effect on the in-phase mode (thin solid line) but 
profoundly affects the out-of-phase mode (thick lines) by 
introducing a gap at w_(fc|| = 0) similar to the neutral 
case. Since in the presence the tunneling the out-of-phase 
mode acquires a gap, the two modes cross at the resonant 
frequency {uj+ — lo-), but the symmetry ("parity") as- 
sociated with the two identical layers does not allow any 
mode coupling or anti-crossing effect. If the two layers in 
the unit cell are not identical then there is a mode cou- 
pling induced anti-crossing around « The inset 
shows the ratio of the oscillator strength of the in-phase 
mode to that of the out-of-phase mode. In sharp contrast 
to the neutral system, in the long wavelength limit the 
out-of-phase mode oj- completely dominates the spectral 
weight in the presence of interlayer tunneling. In the ab- 
sence of tunneling (x = 0), however, the in-phase mode 
0J+ dominates the spectral weight. Our results for the 
collective mode dispersion in the presence of finite single- 
particle tunneling but vanishing Josephson coupling (i.e., 
t ^ 0, Tj — 0) are qualitatively identical to the situation 
with t = 0, Tj 7^ 0, and are therefore not shown sep- 
arately. This is, of course, the expected result because 
t automatically generates an effective Josephson tunnel- 
ing, i.e., an effective Tj, in the superconducting system, 
and therefore the qualitative effect of having a finite Tj 
or a finite t in the superconducting system is similar. 

We also calculate the collective modes of the bilay- 
ered superconducting system by including both the sin- 
gle particle tunneling and the Josephson tunneling be- 
tween the nearest layers (i.e., t,Tj ^ 0). The two lay- 
ers in the unit cell hybridized by the single particle tun- 



neling matrix element, t(k), would lead to a symmet- 
ric and an antisymmetric combination of the quasipar- 
ticlc states for each value of the wave vector k in the 
plane. By introducing the symmetric and antisymmetric 
single electron operators with respect to an interchang- 
ing of the two layers, an,k,cr = -^{cniM + c„2,fco-) and 
l3n,k,a- = -^{cniM - Cn2.ka), the total effective Hamilto- 
nian can be written as 

n kij 

,n'-q + P2, n qT/(q)p2.„'-q] , (27) 



where ei(fc) = e{k) + t{k) and e2(fc) = s{k) — t{k), and 

(28) 

(29) 




and 



^(q) 



u+ U- 
U- u+ 



^(q) 



v+ V- 
V- v+ 



(30) 



where U±^Vi + V2± Tj and V± ^Vi-V2±Tj. This 
Hamiltonian is identical to that in the corresponding two 
subband model, whicb-is well studied in semiconductor 
quantum well systemsEJ. Since there are no off-diagonal 
elements of the interaction with respect to the subband 
index we have well separated intra-subband and inter- 
subband collective modes corresponding to the in-phase 
and out-of-phase modes, respectively. Within our gauge- 
invariant ladder diagram approximation we can easily 
calculate the mode dispersions by following the proce- 
dure outlined in Eqs. (|[) - (p^). The in-phase- mode for 
both the neutral and the charged system is insensitive 
to tunneling in the long wavelength limit, and has essen- 
tially the same long wavelength dispersion as in Eq. ( [2^ ) 
and Eq. ( |25| ) respectively, up to second order in k. The 
out-of-phase mode is, however, strongly affected by both 
the coherent single particle tunneling and the Josephson 
tunneling, and has a dispersion 

Lul {k) ^uji+ [{2tf + vf^k^] [1 + AVo] , (31) 

for neutral superconductors, and 



(1 + Ay^AK)) ^§ + (2^2 + 



UJ_ 



l-^(AF-Al/o) 



(32) 



for charged systems in the presence of finite tunneling. 
In Fig. |3| we show the calculated mode dispersions as a 
function of the in-plane wave vector kud for a fixed 
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FIG. 3. The dispersion of the out-phase mode in the 

charged system in the presence of both the single particle 
tunneling and Josephson tunneling as a function of fcy for a 
fixed fczd = TT. Here, x = Tj/Vo and t is the strength of the 
single particle tunneling with respect to the superconducting 
energy gap. We use the same parameters as Fig. ^ 

kzd = TT. As emphasized before, the collective mode dis- 
persion is quaUtatively independent of the specific tun- 
neling mechanism (i.e., t or Tj), and therefore experi- 
ments involving collective modes cannot distinguish be- 
tween the existing tunneling mechanisms in high- Tc su- 
perconductors as has recently been emphasizedEl in a re- 
lated context. 



III. DISCUSSION AND CONCLUSION 

We calculate in this paper the collective charge density 
fluctuation excitation spectra of both the neutral and the 
charged superconducting bilayerd superlattices with in- 
terlayer intra-cell single particle and Josephson tunnel- 
ing. We use the conserving gauge-invariant ladder dia- 
gram approximation in the Nambu-Gorkov formalism. In 
general, there are two types of density fluctuation modes: 
in-phase {10+) and out-of- phase modes. For neutral 
superconductors, the out-of-phasc collective mode with 
intcrlayer tunneling has a plasma gap depending on the 
tunneling intensity, and the in-phase mode, lying lower 
in energy, dominates the oscillator strength for all wave 
vectors. However, for charged superconductors the two 
phase modes couple to the long range Coulomb interac- 
tion differently, and the out-of-phase mode with tunnel- 
ing dominates the oscillator strength in the long wave- 
length limit — > 0) and finite k^. Since we have used 
two identical 2D layers in each unit cell there is no mode 
coupling effect in our theory between uj± modes at the 
resonant frequency (w-i- ~ lj- ) . If the two layers forming 
the unit cell are not identical then there will be a resonant 



mode coupling effect ("anti-crossing") between the in- 
phase and the out-of-phase modes around w+ « uj- res- 
onance point - the nature of this anti-crossing phenom- 
ena will be similar to what is seen in the corresponding 
intrasubband-intersubband collective—mode coupling in 
semiconductor quantum well systemsEj. We have mostly 
concentrated in the long wavelength regime (fcy ^ 0) - at 
large wave vectors there is significant coupling between 
the collective modes and the pair-breaking excitatior^, 
which has been extensively studied in the literatureBfl. 
We have also neglected the amplitude fluctuation modes 
because they usually carry negligible spectral weights 
compared with the io± phase fluctuation modes. We have 
also used an s-wave ground state symmetry which should 
be a good approximatioio even for d-wave cuprate sys- 
tems as far as the long wavelength collective mode prop- 
erties are concerned. Our use of a BCS-Fermi liquid 
model in our theopv is more difficult to defend except on 
empirical groundsQ and for reasons of simplicity. 

Finally, we consider the effect of intercell tunneling on 
the collective mode spectra, which we have so far ne- 
glected in our consideration. (Our theory includes both 
intracell and intercell Coulomb coupling between all the 
layers, and intracell interlayer single electron and Joseph- 
son tunneling.) The neglect of intercell tunneling is jus- 
tified by the fact that d > c (e.g., in YBCO d = 12 A, 
c — 3A). The general effect of intercell tunneling be- 
comes quite complicated theoretically because one has 
far too many interlayer coupling terms in the Hamilto- 
nian in the presence of both intracell and intercell inter- 
layer tunneling involving both single particle and Joseph- 
son tunneling. It is clear, however, that the main effect 
of a weak intercell interlayer tunneling (either single par- 
ticle or Josephson type, or both) would be to cause a 
2D to 3D transition in the plasma mode by opening up 
a small gap in both uj± modes at long wavelengths (in 
the charged system). The size of this gap (which is the 
effective 3D plasma frequency of the fc^-motion of the 
system) will depend on the intercell tunneling strength. 
This small gap is the 3D c-axis plasma frequency of the 
system, which has bepHpthe subject of several recent stud- 
ies in the literatureoEjiij. 

The introduction of a weak intercell interlayer tunnel- 
ing will therefore modify our calculated results simply 
through a shift of the energy/frequency origin in our cal- 
culated dispersion curves. The origin of the ordinate (i.e., 
the energy/frequency axis) in our results will shift from 
zero to LOc, where lOc is the c-axis plasma frequency aris- 
ing from the intercell interlayer hopping. For an effec- 
tive single band tight binding intercell hopping param- 
eter tc (i.e., the single electron effective bandwidth in 
the c-direction is 2tc), one obtains uJc = coptcd/vp, where 
ujp — [47rnBe^/(Km)]^/^ is the effective 3D plasma fre- 
quency with the 2D a-b plane band mass m [see Eq. 
(p5|)] and vp is the Fermi velocity in the a-b plane. 
Note that uJc ^ tUp because tc is very small by virtue 
of weak intercell coupling. Note also that if one de- 
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fines an efi'ective "3D" c-axis plasma frequency Upc = 
[AimBe^ / {KTric)]^^"^ in analogy with ujp, where rric is now 
the effective mass for electron dynamics along the c-axis, 
then ujc ~ ujpc[t/ {'^■Ef)]^^'^ due to the tight bind nature of 
c-motion. We emphasize that in the presence of intercell 
hopping ojc sets the scale for the lowest energy that a col- 
lective mode can have ii>ittK multilayer superconductoij-^ 
ujc is sometimes referredB'Llu to as a Josephson plasmonO 
in the literature. In general, it is^difficult to theoretically 
estimate ujc in high-Tc materialsO because the effective tc 
(and other parameters) mav-jaflt be known. It is there- 
fore important to emphasizeOO that ujc can be measured 
directly from the c-axis plasma edge in reflectivity ex- 
periments, (we emphasize that a-b plane [plasma edge 
gives ujp and tke c-axis plasma edge gives (jjJIj), and such 
measurementgl3 show that uJc is below the superconduct- 
ing gap in many high- Tc materialsQ. This implies that 
the effective c-axis hopping, tc, in high- Tc materials (ei- 
ther due to single particle hopping or due to Josephson 
coupling arising from coherent Cooper pair hopping) has 
to be very small (much smaller than that given by di- 
rect band structure calculations) in these systems for the 
Josephson plasma frequency lUc to be below the supp.- 
conducting gap, a point first emphasized by Andersonlla. 

The collective mode situation in a bilayer system in 
the presence of both intracell and intercell interlayer cou- 
pling is obviously quite complex, and as emphasized in 
ref. |l^, there could in general be several collective phase 
fluctuation modes depending on the detailed nature of 
intracell and intercell interlayer hopping matrix. In the 
most general bilayer system intercell coupling will give 
rise to two separate w+ plasma bands arising from the 
two distinct possible intercell interlayer coupling — the 
two ujjf- bands lying in energy lower that the two uj— 
bands in the charged systepa. as we show in this paper. In 
the most general situatior£3, there could be two low en- 
ergy Josephson plasma frequencies Wci, ujc2 (>Wci), cor- 
responding to the bottoms of the two uj+ bands, arising 
respectively from the larger and the smaller of the inter- 
cell interlayer hopping amplitudes. To make things really 
complicated one of these modes (cjci) could be below the 
gap and the other (wc2) above the gap, (or, both could be 
below or above the gap). While each of these scenarios 
is possible, c-axis optical response experimental results 
on inter-bilayer charge dynamics in YBCO have been 
interpretedlla to exhibit only one c-axis plasma edge in 
the superconducting state with the frequency lUc between 
60 cm~^ and 200 cm~^, depending on the oxygen con- 
tent. There are three possibilities: (1) The two plasma 
modes (wd ~ Wc2 ~ ^^c) are almost degenerate because 
the corresponding intercell hopping amplitudes are close 
in magnitudes; (2) ujc2 is much lager than cuci (^ ^02) 
because the two intercell hopping amplitude are very dif- 
ferent in magnitudes (we consider this to be an unlikely 
scenario); (3) one of the two modes carries very little 
optical spectral weight and is not showing up in c-axis 
reflectivity measurements, leaving only the other one as 



the observed c-axis plasma edge. There is, in principle, 
a fourth (very unlikely) possibility: the observed plasma 
edge is really llJc2, and the other mode Wci 0Jc2) is too 
low in energy to show up in c-axis reflectivity measure- 
ments. 

Within a nearest-neighbor c-axis interlayer coupling 
model, there is only a single intercell hopping ampli- 
tude, giving rise to only a single c-axis plasma edge lUc, 
which now defines the lowest value that the in-phase col- 
lective mode UJ+ can have, Uc = ujc+ = = 0) — 
LUf. is shifted up from zero at long wavelengths due to 
finite c-axis intercell hopping. The out-of-phase plasma 
edge, Uc- = uj-{k — 0), will obviously lie much higher 
in energy than lUc+ = lUc because the intracell inter- 
layer hopping is much stronger than the intercell in- 
terlayer hopping. In particular, even t.hnngh the lUc+ 
mode may lie in the superconducting gaplljllj, we expect 
LUc- to lie much above the superconducting gap energy 
in YBCO. A crude qualitative estimate can be made 
by assuming that the intra- and intercell hopping am- 
plitudes scale as inverse squares of lattice parameters: 
^intra/^intor ~ (d/c)'^ = 16. This then Icads to the ap- 
proximate formula Wc- ~ 16^ LiJc+ = 256 ujc, which, for 
YBCO, implies that the long wavelength out-of-phase 
mode should lie between 2 eV and 6 eV, depeding on the 
oxygen content (assuming that the c-axis plasma edge 
varies between 60 cm^^ and 200 cm^^, as reported in 
ref. |l^, depending on the oxygen content). While there 
is some minor observable structure in optical experiments 
at high energies, we cannot find any compelling evidence 
in favor of the existence of a high energy out-of-phase 
mode in the currently available experimental data. We 
feel that a spectroscopic experiment, using, for exam- 
ple, the inelastic electron energy loss spectroscopy which 
could probe the mode dispersion (and which has been 
highly successful in studying bulk plasmons in metal 
films) of the uj- mode at high energy, may be required to 
unambiguously observe the out-of-phase collective mode. 
What we have shown in this paper is that under suitable 
conditions (finite k and fc^) the out-of-phase mode 
carries reasonable spectral weight and should be observ- 
able in principle — actual observation, however, awaits 
experimental investigations using external probes which 
can study mode dispersion at finite wave vectors (which 
optical experiments by definition cannot do; they are long 
wavelength probes). 
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